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We review various issues related to the direct detection of constituents of 
dark matter, which are assumed to be Weakly Interacting Massive Particles 
(WIMPs). We specifically consider heavy WIMPs such as: 1) The lightest su- 
persymmetric particle LSP or neutralino. 2) The lightest Kaluza-Klein particles 
in theories of extra dimensions and 3) other extensions of the standard model. 
In order to get the event rates one needs information about the structure of 
the nucleon as well as as the structure of the nucleus and the WIMP velocity 
distribution. These are also examined Since the expected event rates for de- 
tecting the recoiling nucleus are extremely low and the signal does not have a 
characteristic signature to discriminate them against background we consider 
some additional aspects of the WIMP nucleus interaction, such as the periodic 
behavior of the rates due to the motion of Earth (modulation effect). Since, 
unfortunately, this is characterized by a small amplitude we consider other op- 
tions such as directional experiments, which measure not only the energy of 
the recoiling nuclei but their direction as well. In these, albeit hard, experi- 
ments one can exploit two very characteristic signatures: a)largc asymmetries 
and b) interesting modulation patterns. Furthermore we extended our study to 
include evaluation of the rates for other than recoil searches such as: i) Transi- 
tions to excited states, ii) Detection of recoiling electrons produced during the 
ncutralino-nucleus interaction and iii) Observation of hard X-rays following the 
de-excitation of the ionized atom. 

Keywords: WIMP; Dark Matter; CDM; Neutralino; Supersymmetry; LSP; 
Kaluza-Klein WIMPs; Modulation; Directional event rate. 

1. Introduction 

The combined MAXIMA-1 1 , BOOMERANG 2 , DASI 3 , COBE/DMR Cos- 
mic Microwave Background (CMB) observations 4 , the recent WMAP data 5 
and SDSS 6 imply that the Universe is flat 7 and that most of the matter in 
the Universe is dark, i.e. exotic. These results have been confirmed and im- 
proved by the recent WMAP data. 8 The deduced cosmological expansion is 
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consistent with the luminosity distance as a function of redshift of distant 
supernovae. 9-11 According to the scenario favored by the observations there 
are various contributions to the energy content of our Universe. The most 
accessible energy component is baryonic matter, which accounts for ~ 5% 
of the total energy density. A component that has not been directly ob- 
served is cold dark matter (CDM)): a pressureless fluid that is responsible 
for the growth of cosmological perturbations through gravitational insta- 
bility. Its contribution to the total energy density is estimated at ~ 25%. 
The dark matter is expected to become more abundant in extensive halos, 
that stretch up to 100-200 kpc from the center of galaxies. The component 
with the biggest contribution to the energy density has an equation of state 
similar to that of a cosmological constant and is characterized as dark en- 
ergy. The ratio w = p/p is negative and close to — 1. This component is 
responsible for ~ 70% of the total energy density and induces the observed 
acceleration of the Universe 9-11 . The total energy density of our Universe 
is believed to take the critical value consistent with spatial flatness. Addi- 
tional indirect information about the existence of dark matter comes from 
the rotational curves 12 . The rotational velocity of an object increases so 
long is surrounded by matter. Once outside matter the velocity of rotation 
drops as the square root of the distance. Such observations are not possi- 
ble in our own galaxy. The observations of other galaxies, similar to our 
own, indicate that the rotational velocities of objects outside the luminous 
matter do not drop. So there must be a halo of dark matter out there. 

Since the non exotic component cannot exceed 40% of the CDM 13 , 
there is room for exotic WIMP's (Weakly Interacting Massive Particles). 
In fact the DAMA experiment 14 has claimed the observation of one signal in 
direct detection of a WIMP, which with better statistics has subsequently 
been interpreted as a modulation signal 15 . These data, however, if they 
are due to the coherent process, are not consistent with other recent exper- 
iments, see e.g. EDELWEISS and CDMS 16 . It could still be interpreted as 
due to the spin cross section, but with a new interpretation of the extracted 
nucleon cross section. 

Since the WIMP is expected to be very massive, m x > 30GeU, and 
extremely non relativistic with average kinetic energy T < lOOKeV, it can 
be directly detected mainly via the recoiling of a nucleus in the WIMP- 
nucleus elastic scattering. 

The above developments are in line with particle physics considerations. 

(1) Dark matter in supersymmetric theories 

The lightest supersymmetric particle (LSP) or ncutralino is the most 
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natural WIMP candidate. In the most favored scenarios the LSP can 
be simply described as a Majorana fermion, a linear combination of the 
neutral components of the gauginos and Higgsinos 12-17 . 
In order to compute the event rate one needs an effective Lagrangian 
at the elementary particle (quark) level obtained in the framework of 
supersymmetry 12 > 17 > 18 . One starts with representative input in the 
restricted SUSY parameter space as described in the literature, e.g. 
Ellis et al 19 , Bottino et al , Kane et al , Castano et al and Arnowitt 
et al 18 as well as elsewhere 20 ~ 21 . We will not, however, elaborate on 
how one gets the needed parameters from supersymmetry. 
Even though the SUSY WIMPs have been well studied, for tor the 
reader's convenience we will give a description in sec. 2 of the basic 
SUSY ingredients needed to calculate LSP-nucleus scattering cross sec- 
tion. 

(2) Kaluza-Klcin (K-K) WIMPs. 

These arise in extensions of the standard model with compact extra di- 
mensions. In such models a tower of massive particles appear as Kaluza- 
Klein excitations. In this scheme the ordinary particles are associated 
with the zero modes and are assigned K-K parity +1. In models with 
Universal Extra Dimensions one can have cosmologically stable parti- 
cles in the excited modes because of a discreet symmetry yielding K-K 
parity —1 (see previous work 22-24 as well as the recent review by Ser- 
vant 25 ). 

The kinematics involved is similar to that of the neutralino, leading to 
cross sections which are proportional /z 2 , \i r being the WIMP-nucleus 
reduced mass. Furthermore the nuclear physics input is independent of 
the WIMP mass, since for heavy WIMP mu r ~ Am p . There are appear 
two differences compared to the neutralino, though, both related to its 
larger mass. 

i) First the density (number of particles per unit volume) of a WIMP 
falls inversely proportional to its mass. Thus, if the WIMP's considered 
are much heavier than the nuclear targets, the corresponding event rate 
takes the form: 

R(m WIMP ) = R(A) A GCV (1) 
mwiMP 

where R(A) are the rates extracted from experiment up to WIMP 
masses of the order of the mass of the target. 

ii) Second the average WIMP energy is now higher. In fact one finds 
that (Twimp) — jmwiMpVQ ~ 40(m W / M p/(100 GeV))keV (v ~ 
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2.2 x 10 5 km/s). Thus for a K-K WIMP with mass 1 TeV, the average 
WIMP energy is 0.4 McV. Hence, due to the high velocity tail of the 
velocity distribution, one expects an energy transfer to the nucleus 
in the MeV region. Thus many nuclear targets can now be ex- 
cited by the WIMP-nucleus interaction and the de-excitation 
photons can be detected. 

In addition to the particle model one needs the following ingredients: 

• A procedure in going from the quark to the nucleon level, i.e. a quark 
model for the nucleon. The results depend crucially on the content of 
the nucleon in quarks other than u and d. This is particularly true for 
the scalar couplings as well as the isoscalar axial coupling 26-27 . Such 
topics will be discussed in sec. 4. 

• computation of the relevant nuclear matrix elements 28 ~~ 29 using as re- 
liable as possible many body nuclear wave functions. By putting as 
accurate nuclear physics input as possible, one will be able to constrain 
the SUSY parameters as much as possible. The situation is a bit sim- 
pler in the case of the scalar coupling, in which case one only needs the 
nuclear form factor. 

• Convolution with the LSP velocity Distribution. To this end we will 
consider here Maxwell-Boltzmann 12 (MB) velocity distributions, with 
an upper velocity cut off put in by hand. The characteristic velocity of 
the M-B distribution can be increased by a factor n (vq — > hvq, n > 
l)by considering the interaction of dark matter and dark energy. 30 
Other distributions are possible, such as non symmetric ones, like those 
of Drukier 31 and Green 32 , or non isothermal ones, e.g. those arising 
from late in-fall of dark matter into our galaxy, like Sikivie's caustic 
rings 33 . In any event in a proper treatment the velocity distribution 
ought to be consistent with the dark matter density as, e.g., in the 
context of the Eddington theory 34 . 

Since the expected rates are extremely low or even undetectable with 
present techniques, one would like to exploit the characteristic signatures 
provided by the reaction. Such are: 

(1) The modulation effect, i.e the dependence of the event rate on the 
velocity of the Earth 

(2) The directional event rate, which depends on the velocity of the sun 
around the galaxy as well as the the velocity of the Earth, has recently 
begun to appear feasible by the planned experiments 35,36 . 
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(3) Detection of signals other than nuclear recoils, such as 

• Detection of 7 rays following nuclear de-excitation, whenever 
possible 37,38 . This seems to become feasible for heavy WIMPs 
especially in connection with modified M-B distributions due 
to the coupling of dark matter and dark energy ((Twimp) — 
n 2 40 (m WIM p/(l00 GeV)) , n > IkeV) 

• Detection of ionization electrons produced directly in the LSP- 
nucleus collisions 39,40 . 

• Observations of hard X-rays produced 41 , when the inner shell 
electron holes produced as above are filled. 

In all calculations we will, of course, include an appropriate nuclear 
form factor and take into account the influence on the rates of the detector 
energy cut off. We will present our results a function of the LSP mass, m x , 
in a way which can be easily understood by the experimentalists. 

2. The Feynman Diagrams Entering the Direct Detection 
of WIMPS. 

2.1. The Feynman Diagrams involving the neutralino 

The neutralino is perhaps the most viable WIMP candidate and has been 
extensively studied (see, e.g., our recent review 42 ). Here we will give a very 
brief summary of the most important aspects entering the direct neutralino 
searches. In currently favorable supergravity models the LSP is a linear 
combination 12 of the neutral four fermions B, W3, Hi and Hi which are the 
supersymmetric partners of the gauge bosons and and the Higgs 
scalars Hi and i?2- Admixtures of s- neutrinos are expected to be negligible. 
The relevant Feynman diagrams involve Z-exchange, s-quark exchange and 
Higgs exchange. 

2.1.1. The Z-exchange contribution. 

The relevant Feynman diagram is shown in Fig. 1. It does not lead to 
coherence, since ^7^^ = for a Majorana fermion like the neutralino (the 
Majorana fermions do not have electromagnetic properties). The coupling 
^7a75* yields negligible contribution for a non relativistic particle in the 
case of the spin independent cross section. 43 It may be important in the 
case of the spin contribution, which arises from the axial current). 
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Fig. 1. The LSP-quark interaction mediated by Z and Higgs exchange. 
2.1.2. The s-quark Mediated Interaction 

The other interesting possibility arises from the other two components of 
Xi, namely B and W3. Their corresponding couplings to s-quarks (see Fig. 
2 ) can be read from the appendix C4 of Ref. 17 and our earlier review. 42 
Normally this contribution yields vector like contribution, i.e it does not 




Fig. 2. The LSP-quark interaction mediated by s-quark exchange. Normally it yields 
V-A interaction which does not lead to coherence at the nuclear level. If, however, the 
isodoublet s-quark is admixed with isosinglet one to yield a scalar interaction at the 
quark level. 

lead to coherence. If, however, there exists mixing between the s-quarks 
with isospin 1/2 (qx) and the isospin (qn), the s-quark exchange may 
lead to a scalar interaction at the quark level and hence to coherence over 
all nucleons at the nuclear level. 42 

2.1.3. The Intermediate Higgs Contribution 

The most important contribution to coherent scattering can be achieved 
via the intermediate Higgs particles which survive as physical particles. In 
supersymmetry there exist two such physical Higgs particles, one light h 
with a mass <120 GeV and one heavy H with mass m#, which is much 
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larger. The relevant interaction can arise out of the Higgs-Higgsino-gaugino 
interaction 42 leading to a Feynman diagram shown in Fig. 1. 

In the case of the scalar interaction the resulting amplitude is propor- 
tional to the quark mass. 

2.2. The Feynman Diagrams involving the K-K WIMPs 

2.2.1. The Kaluza-Klein Boson as a dark matter candidate 

We will assume that the lightest exotic particle, which can serve as a dark 
matter candidate, is a gauge boson B 1 having the same quantum numbers 
and couplings with the standard model gauge boson B, except that it has 
K-K parity —1. Thus its couplings must involve another negative K-K parity 
particle. In this work we will assume that such a particle can be one of the 
K-K quarks, partners of the ordinary quarks, but much heavier 22-2 . 

• Intermediate K-K quarks. 

this case the relevant Feynman diagrams are shown in fig. 3. 



q B m 



Ji) 



BO) 



Fig. 3. K-K quarks mediating the interaction of K-K gauge boson B with quarks at 
tree level. 



The amplitude at the nucleon level can be written as 
Meoh = A(e*'.e)N 



11 + 12t 3 m p m w 1 + t 3 m w 

h(A-) "I o M^J 



54 (m B( i)) 5 



3 m B ( 



1 + A + A 2 12 

A = ikJ^Gprnw tan 2 6> w ,/i(A) = A2 ( 1 + A / 2 ) 2 . 



N 
(2) 



/a (A) 



1 + A A _ m q( i } _ i 



A(l + A/2) 



We see that the amplitude is very sensitive to the parameter A ("res- 
onance effect"). 
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In going from the quark to the nucleon level the best procedure is 
to replace the quark energy by the constituent quark mass ~ l/3m p , 
as opposed to adopting 22-24 a procedure related to the current mass 
encountered in the neutralino case. 42 

In the case of the spin contribution we find at the nucleon level that: 

M spin = -UV2Gf toh/ tan 2 9 W \ m P mw _ / 1 (A)j(e*' x e). 

[N<r(g + gi T 3 )N] (3) 



17 A A 7 A 17 A A 7 

on = — Am H Ad H As , q 1 = — Au Ad 

wu 18 18 18 ' y 18 18 

for the isoscalar and isovector quantities. 42 The quantities A q are given 
by 42 

Aw = 0.78 ± 0.02 , Ad = -0.48 ± 0.02 , As = -0.15 ± 0.02 

We thus find g = 0.26 , g x = 0.41 => a p = 0.67 , a n = -0.15. 
• Intermediate Higgs Scalars. 

The corresponding Feynman diagram is shown in Fig. 4 The relevant 




Fig. 4. The Higgs H mediating interaction of K-K gauge boson B 1 with quarks at tree 
level (on the left). The Z-boson mediating the interaction of K-K neutrino i/W with 
quarks at tree level (on the right). 



amplitude is given by: 

M N (h) = -i 4v / 2G F TO^tan 2 6» H / 



(4) 

In going from the quark to the nucleon level we follow a procedure 
analogous to that of the of the neutralino, i.e. -< N\m q qq\N f q m p 
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2.2.2. K-K neutrinos as dark matter candidates 

The other possibility is the dark matter candidate to be a heavy K-K neu- 
trino. We will distinguish the following cases: 

• Process mediated by Z-exchange. 

The amplitude associated with the diagram of Fig. 4 becomes: 

M vm = -^=G F J\^)J X (NNZ) (5) 

with J\(NNZ) the standard nucleon neutral current and 

Ja(^ (1) ) = ^ (1) 7a75^ (1) , Ja(^ (1) )=^ (1) 7a(1-75)^ (1) 

for Majorana and Dirac neutrinos respectively. 

• Process mediated by right handed currents via Z'-boson exchange. 
The process is similar to that exhibited by Fig. 4, except that instead of 
Z we encounter Z', which is much heavier. Assuming that the couplings 
of the Z' are similar to those of Z, the above results apply except 
that now the amplitudes are retarded by the multiplicative factor k = 

2 / 2 

m z /m z , 

• Process mediated by Higgs exchange. 

In this case in Fig 4 the Z is replaced by the Higgs particle. Proceeding 
as above we find that the amplitude at the nucleon level is: 

M vW (h) = -2V2G F ^^M V W <N\N^Y] f q (6) 
m h q 

In the evaluation of the parameters f q one encounters both theoretical 
and experimental errors. 

3. Other non SUSY Models 

We should mention that there exist extensions of the standard model not 
motivated by symmetry. Such are: 

• Models which introduce extra higgs particles and impose a discrete 
symmetry which leads to a "parity" a la R-parity or K-K parity. 44 

• Extensions of the standard model, which do not require a parity, but 
introduce high weak isospin multiplcts 45 with Y=0. So the WIMP- 
nucleus interaction via Z-exchange at tree level is absent and the dom- 
inant contribution to the WIMP-nucleus scattering occurs at the one 
loop level. 
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• Another interesting extension of the standard model is in the direction 
of tecnicolor. 46 In this case the WIMP is the neutral LTP (lightest neu- 
tral technibaryon) . This is scalar particle, which couples to the quarks 
via derivative coupling through Z-exchange. 

4. Going from the Quark to the Nucleon Level 

In going from the quark to the nucleon level one has to be a bit more careful 
in handling the quarks other than u and d. This is especially true in the 
case of the scalar interaction, since in this case the coupling of the WIMP 
to the quarks is proportional to their mass 42 . Thus one has to consider in 
the nucleon not only sea quarks (uu, dd and ss) but the heavier quarks as 
well due to QCD effects 47 . This way one obtains the scalar Higgs-nucleon 
coupling by using effective parameters f q defined as follows: 



where is the nucleon mass. The parameters f g , q = u,d,s can be 
obtained by chiral symmetry breaking terms in relation to phase shift and 
dispersion analysis (for a recent review see 42 ). We like to emphasize here 
that since the current masses of the u and d quarks are small, the heavier 
quarks tend to dominate even though the probability of finding them in 
the nucleus is quite small. In fact the s quark contribution may become 
dominant, e.g. allowed by the above analysis is the choice: 



f d = 0.046, /„ = 0.025, f s = 0.400, f c = 0.050, f b = 0.055, f t = 0.095 



The isoscalar and the isovector axial current in the case of K-K theo- 
ries has already been discussed above. In the case of the neutralino these 
couplings at the nucleon level, f\, are obtained from the correspond- 
ing ones given by the SUSY models at the quark level, /^(<z), f\ (q), via 
renormalization coefficients g% g\, i.e. f\ = S^/aO?)* /A = fli/AO?)- Thc 
renormalization coefficients are given terms of Aq defined above, 48 via the 
relations 

g° A = Au + Ad + As = 0.77 - 0.49 - 0.15 = 0.13 , g\ = Au - Ad = 1.26 

We see that, barring very unusual circumstances at thc quark level, the 
isoscalar contribution is negligible. It is for this reason that one might prefer 
to work in the isospin basis. 




(7) 
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5. The allowed SUSY Parameter Space 

It is clear from the above discussion that the LSP-nucleon cross section 
depends, among other things, on the parameters of supcrsymmctry. One 
starts with a set of parameters at the GUT scale and predicts the low energy 
observables via the renormalization group equations (RGE). Conversely 
starting from the low energy phenomenology one can constrain the input 
parameters at the GUT scale. The parameter space is the most crucial. In 
SUSY models derived from minimal SUGRA the allowed parameter space 
is characterized at the GUT scale by five parameters: 

• two universal mass parameters, one for the scalars, to , and one for the 
fermions, m^. 

• tan/3, i.e the ratio of the Higgs expectation values, (H 2 ) / (Hi). 

• The trilinear coupling A (or m\ ole ) and 

• The sign of n in the Higgs self-coupling uHiH 2 - 

The experimental constraints 42 restrict the values of the above parameters 
yielding the allowed SUSY parameter space. 



6. Event rates 

The differential non directional rate can be written as 

p(0) m , . ... 
dRundir = -; da(u,v)\v\ (8) 

where A is the nuclear mass number, p(0) w 0.3GeV/cm 3 is the WIMP 
density in our vicinity, m is the detector mass, m x is the WIMP mass and 
da(u,v) is the differential cross section. 

The directional differential rate, i.e. that obtained, if nuclei recoiling in 
the direction e are observed, is given by 42 ' 49 : 

p(0) m „ 1 \A* ... 
dRdir = -. \v\v.e Bfw.e) — da(u,v of = — v.e) 

where 0(x) is the Heaviside function. 

The differential cross section is given by: 

du — — 
da( U , V) == — [(S S F(U) 2 + Eapin^ 11 («)] (9) 

where u the energy transfer Q in dimensionless units given by 

U = 7T > Qo = [m p Ab]- 2 = A0A- A/3 MeV (10) 
Qo 
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with b is the nuclear (harmonic oscillator) size parameter. F(u) is the nu- 
clear form factor and Fn(u) is the spin response function associated with 
the isovector channel. 

The scalar contribution is given by: 

2 



2„S' 



1 



fs 2Z-A 
7j A 

_1_ Ik. 



Hr(p) 



) 2 A 2 



(11) 



(since the heavy quarks dominate the isovector contribution is negligible). 



is the LSP-nucleon scalar cross section. 



The spin contribution is given by: 

^ s rn I 



Mr \2 spin > > 

' ^v.\° ^spini^spin 



>r(p) 



3(1 + 4) 2 

•> A 



S{u) 



(12) 



S{u) 



5(0) = [An (o)) 2 + 2^fi (o)fii(o) + ^(o)) 2 ; 

J A 



(13) 



The couplings f\ (f\) and the nuclear matrix elements fii (0) (f2o(0)) as- 
sociated with the isovector (isoscalar) components are normalized so that, 
in the case of the proton at u = 0, they yield ( S pin = 1- 

With these definitions in the proton neutron representation we get: 



{spin = -S (0) , S (0) = 



(^„(0)) 2 + 2^^(0)^(0) + ft 2 (0) 



(14) 



where 0^,(0) and Sl„ (0) are the proton and neutron components of the 
static spin nuclear matrix elements. In extracting limits on the nuclcon 
cross sections from the data we will find it convenient to write: 



spin > 
*^P>X° ^ s P* n 



j—— ^ ^rj 



io n e 



iS\2 



(15) 



In Eq. (15) S the relative phase between the two amplitudes a p and a n , 
which in most models is or n, i.e. one expects them to be relatively real. 
The static spin matrix elements are obtained in the context of a given 
nuclear model. Some such matrix elements of interest to the planned ex- 
periments can be found in. 42 

The spin ME arc defined as follows: 



n p (o) 



J+l 



-< JJ\a z (p)\J Jy , 0,(0) 



J + l 



-< J J\a z {ri)\J J y 
(16) 
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where J is the total angular momentum of the nucleus and o z = 2S Z . 
The spin operator is defined by S z (p) — J2i=i Sz{i), i-e. a sum over all 
protons in the nucleus, and S z (n) — J2iLi ^z(i), i- e - a sum over all neutrons. 
Furthermore fi (0) = H P (0) + fi„(0) , fii(0) = fi p (0) - fi n (0) 



7. The WIMP velocity distribution 

To obtain the total rates one must fold with WIMP velocity distribution 
and integrate the above expressions over the energy transfer from Q m in 
determined by the detector energy cutoff to Q max determined by the max- 
imum LSP velocity (escape velocity, put in by hand in the Maxwellian 
distribution), i.e. v esc — 2.84 w with v the velocity of the sun around the 
center of the galaxy (229 Km/ s). 

For a given velocity distribution f(u'), with respect to the center of the 
galaxy, one can find the velocity distribution in the Lab {(v,ve) by writing 
v = v + ve , ve=v + vi, with v i the Earth's velocity around the sun. 

It is convenient to choose a coordinate system so that x is radially out 
in the plane of the galaxy, z in the sun's direction of motion and y = z x ob. 

Since the axis of the ecliptic lies very close to the x, y plane (u) — 186.3°) 
only the angle 7 = 29.8° becomes relevant. Thus the velocity of the earth 
around the sun is given by 

ve = Vqz + V\ ( since x — cosa cos-f y + cosa siwy z ) (17) 

where a is phase of the earth's orbital motion. 

The WIMP velocity distribution i{v') is not known. Many velocity dis- 
tributions have been used. The most common one is the M-B distribution 
with characteristic velocity Vq with an upper bound v esc = 2.84u - 

f( v >) = - e -(V/v ) 2 (18) 

Modifications of this velocity distribution have also been considered such 
as: i) Axially symmetric M-B distribution. 31 ' 50 and ii) modifications of the 
characteristic parameters of the M-B distribution by considering a coupling 
between dark matter and dark energy 30 (vo — > nvo,v esc nv esc ). Other 
possibilities are adiabatic velocity distribution following the Eddington ap- 
proach 51-52 , caustic rings 53-54 and Sagittarius dark matter 32 . 
For a given energy transfer the velocity v is constrained to be 



QAm p 1 
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8. The Direct detection rate 

The event rate for the coherent WIMP-nuclcus elastic scattering is given 

by. 42, 49, 55, 56 

p(0) m 



R 
with 



VW) fcoh(A,Hr(A))*2 iX o +f Sp in(A, t l r (A))a S p P ™ (. 



^sptn 



(20) 



fcoh(A, Hr{A)) 



lOOGeV 



fj, r (A) 



1 2 



A t coh (1 + h coh cosa) (21) 



fspin(A, /J, r (A)) 



Mr(p) 



,(A) 



t sp in (1 + hspincosa) (22) 



with cr^ x „ and f* 3 ^" the scalar and spin proton cross sections (spin the 
nuclear spin ME. In the above expressions h is the modulation amplitude. 

The number of events in time t due to the scalar interaction, which leads 
to coherence, is: 



R - 1 - 6010 l7o.3G e y C m-3T^280km S -i 10^ -g/cofc^W^)) 

(23) 

In the above expression m is the target mass, A is the number of nucleons 
in the nucleus and (v 2 ) is the average value of the square of the WIMP 
velocity. 

In the case of the spin interaction we write: 



i?~ 16: 



P(0) 



a p,x° 



ly0.3GeVcm-* lK g 280km S -i 10~ 2 p b-/WAM^)) (24) 

Note the different scale for the proton spin cross section. The parameters 
fcoh{A, n r (A)), f sp in(A, fi r (A)), which give the relative merit for the co- 
herent and the spin contributions in the case of a nuclear target compared 
to those of the proton, have already been tabulated 42 for energy cutoff 
Qmin = 0, 10 keV. It is clear that for large A the coherent process is ex- 
pected to dominate unless for some reason the scalar proton cross section 
is very suppressed. 

In the case of directional experiments the event rate is given by Eqs (23) 
and (24) except that now: 

2 



fcoh{A,Hr(A)) 



lOOGeV 



Hr (A) 



A—t coh (1 + h m (coh)cos(a + a m n)) 
Ztt 

(25) 
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(1 + h m (spin)cos(a + a m ir)) 
(26) 

In the above expressions h m is the modulation amplitude and a m the shift 
in the phase of the modulation (in units of it) relative to the phase of 
the Earth. k/(2tt), k < 1, is the suppression factor entering due to the 
restriction of the phase space, k, h m and a m depend on the direction of 
observation. It is precisely this dependence as well as the large values of h m , 
which can be exploited to reject background, 42 that makes the directional 
experiments quite attractive in spite of the suppression factor relative to 
the standard experiments. 



lOOGcV fi r (A) 



9. Bounds on the scalar proton cross section 

Using the above formalism one can obtain the quantities of interest t and 
h both for the standard as well as the directional experiments. Due to 
lack of space we are not going to present the obtained results here. The 
interested reader can find some of these results elsewhere 42,49 . Here we 
are simply going to show how one can employ such results to extract the 
nucleon cross section from the data. Due to space considerations we are 
not going to discuss the limits extracted from the data on the spin cross 
sections, since in this case one has to deal with two amplitudes (one for the 
proton and one for the neutron). We will only extract some limits imposed 
on the scalar nucleon cross section (the proton and neutron cross section 
are essentially the same). In what follows we will employ for all targets 57-58 
the limit of CDMS II for the Ge target 59 , i.e. < 2.3 events for an exposure 
of 52.5 Kg-d with a threshold of 10 keV. This event rate is similar to that for 
other systems. 60 The thus obtained limits are exhibited in Fig. 5. For larger 
WIMP masses one can extrapolate these curves, assuming an increase as 

10. Transitions to excited states 

The above formalism can easily be extended to cover transitions to excited 
states. Only the kinematics and the nuclear physics is different. In other 
words one now needs: 

• The inelastic scalar form factor. 

The transition amplitude is non zero due to the momentum transfer 
involved. The relevant multipolarities are determined by the spin and 
parity of the final state. 
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Fig. 5. The limits on the scalar proton cross section for A= 127 on the left and A= 73 
on the right as functions of m x . The continuous (dashed) curves correspond to Q m in = 
(10) keV respectively. Note that the advantage of the larger nuclear mass number of 
the A= 127 system is counterbalanced by the favorable form factor dependence of the 
A= 73 system. 



• Spin induced transitions. 

In this case one can even have a Gamow- Teller like transition, if the 
final state is judiciously chosen. 

In the case of 127 1 the static spin matrix element involving the first 
excited state around 50 keV is twice as large compared to that of the 
ground state. 38 The spin response function was assumed to be the same 
with that of the ground state. The results obtained 38 are shown in Fig. 
6. These results are very encouraging, since, as we have mentioned, for 




m L sp — > (GeV) 



Fig. 6. The ratio of the rate to the excited state divided by that of the ground state as a 
function of the LSP mass (in GeV) for 127 1. We found that the static spin matrix element 
of the transition from the ground to the excited state is a factor of 1.9 larger than that 
involving the ground state and assumed that the spin response functions Fn(u) are the 
same. On the left we show the results for Q m i n = and on the right for Q m i n = 10 KeV. 



heavier WIMPS like those involved in K-K theories, the branching ratios 
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are expected to be much larger. Thus one may consider such transitions, 
since the detection of de-excitation 7 rays is much easier than the detection 
of recoiling nuclei. 

11. Other non recoil experiments 

As we have already mentioned the nucleon recoil experiments are very hard. 
It is therefore necessary to consider other possibilities. One such possibility 
is to detect the electrons produced during the WIMP-nucleus collisions 39 ' 40 
employing detectors with low energy threshold with a high Z target. Better 
yet one may attempt to detect the very hard X-rays generated when the 
inner shell electron holes arc filled. The relative X-ray to nucleon recoil 
probabilities [Za K /a r ]i, for i = L(m x < lOOGeV), M(100 GeV < m x < 
200 GeV) and H(m x ~ 200 GeV) are shown in table 1. For even heavier 
WIMPs, like those expected in K-K theories, the relative probability is 
expected to be even larger. 



K X-ray 


E K (K t3 ) keV 




,Z<T K {K i:j ), 

1 (7 r \m 




K a2 


29.5 


0.0086 


0.0560 


0.0645 


K a i 


29.8 


0.0160 


0.1036 


0.1196 


K/31 


33.6 


0.0047 


0.0303 


0.0350 


K fj2 


34.4 


0.0010 


0.0067 


0.0077 



The K Q and lines can be separated experimentally by using good 
energy-resolution detectors, but the sum of all K lines can be measured in 
modest energy-resolution experiments. 

12. Conclusions 

We examined the various signatures expected in the direct detection of 
WIMPs via their interaction with nuclei. We specially considered WIMPs 
predicted in supersymmetric models (LSP or neutralino) as well as theo- 
ries with extra dimensions. We presented the formalism for the modulation 
amplitude for non directional as well as directional experiments. We dis- 
cussed the role played by nuclear physics on the extraction of the nucleon 
cross sections from the data. We also considered non recoil experiments, 
such as measuring the 7 rays following the de-excitation of the nucleus 
and/or the hard X-rays after the de-excitation of the inner shell electron 
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holes produced during the WIMP nucleus interaction. These are favored by 
very heavy MIMPs in the TeV region and velocity distributions expected 
in models allowing interaction of dark matter and dark energy. 
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